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Abstract 

The linear delta expansion is applied to a calculation of the SU{2) mass gap on 
the lattice. Our results compare favourably with the strong-coupling expansion and 
are in good agreement with recent Monte Carlo estimates. 



1. Introduction 

The linear delta expansion (LDE) is an analytic approach to field theory which has been 
applied to a number of different problems (see for example Ref. The approach is non- 
perturbative in the sense that a power series expansion is made in a parameter 5 artificially 
inserted into the action, rather than in a coupling constant of the theory. The calculational 
techniques required do not differ greatly from conventional Feynman diagrams. An essential 
part of the approach is an optimization with respect to another parameter, in the present 
case J, appearing in the ^-extended action. 

The linear delta expansion uses 5 as an interpolation between a soluble action So and 
the action for the desired theory S. The action is written: 

Ss = (l-5)S + 5S, (1) 

where So contains some dependence on the optimization parameter J. A vacuum generating 
functional or appropriate Green function may then be evaluated as a power series in S, 
which is set equal to unity at the end of the calculation. 

Of course this power series is only calculated to a finite number of terms, and will 
therefore retain some dependence on J which would be absent in the sum to all orders 
when 6 is set equal to one. A well-motivated criterion for fixing J is to demand that, 
at least locally, the truncated result should be independent of J. This is the principle of 
minimal sensitivity (PMS) ||. If Cat denotes the iVth approximant to a quantity C, the 
requirement is 

This, or some similar criterion, is an intrinsic part of the LDE, providing the non- 
perturbative dependence on the coupling constant of the theory. For example, in the delta 
expansion of the integral / <ia;exp(— gx 4 ), the PMS correctly reproduces its g~ x ^ depen- 
dence. The application of the PMS is also vital for the convergence of the 5 series, which 



has been rigorously proved for the zero- dimensional 4 vacuum generating functional || 
and the finite temperature partition function of the anharmonic oscillator in quantum me- 
chanics pi. The proof has been recently extended H to the connected vacuum generating 
function W = In Z in zero dimensions. 

A number of non-perturbative approaches to field theory are related to the LDE. At 
first order in 6, the LDE is related to the Gaussian approximation ||, and at higher orders 
to generalizations of this |7j. It also has much in common with work of Kleinert || and of 
Sissakian et al. ||. 

In the context of lattice gauge theories the LDE has been applied, with various choices 



of the trial action |10]-(T2|], to the groups Z%, U(l) and SU(2), mainly in calculating the 
plaquette energy E P . A particularly useful trial action is the one proposed by Zheng et 
al. [13[], based on single links: 

S = Jj2^U e (3) 

i 

The sum runs over all links t of the lattice, and the parameter J is used for optimization. 
However, these authors used a different optimization criterion, more closely related to the 
conventional variational method, in which a rigorous inequality for the free energy at 0(5) 
was applied at all orders in 5. Such a procedure is liable to forfeit the convergence which 
may be provided by an order-by-order optimization. 

Two of the present authors |T3| used the Zheng trial action with the PMS in its usual 



sense in a calculation of the SU{2) plaquette energy. This was found to give excellent 
agreement at 0(5 3 ) with Monte Carlo results in the weak coupling regime. Following on 
from this, the phase structure of the mixed SU{2) - SO (3) action was studied |15| , and 
again the results to 0(5 3 ) gave good agreement with the Monte Carlo results. 

We were therefore encouraged to attempt to extend the method to the more difficult 
problem of the mass gap in lattice SU (2) using the same trial action. For such quantities, 
which involve finding the exponential fall-off of a correlator at large separations, semi- 



analytic methods have, in principle, an advantage over Monte-Carlo methods, insofar as 
the size of the lattice is not limited and small signals are not masked by statistics. 

In Section 2 we set up the formalism for the problem to be studied and explain how 
the diagrams which arise in the delta expansion of the modified action are evaluated. The 
optimization procedure adopted is explained in Section 3, where the results are presented 
first in lattice units, and then in terms of the SU(2) lattice constant by looking for the 
correct scaling limit as a — > 0. In Section 4 we summarize the paper and indicate some 
directions for further development. The appendix shows how the evaluation of expectation 
values in the background of the trial action of Eq. (Q) can be organized in a way amenable 
to symbolic computation. 
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2. The Diagrammatic Expansion 

We consider the SU (2) gauge theory on the lattice. The ^-extended action is : 

S& = Ufi^Up + (1 - 6)J^2tTU t (4) 
^ p t 

The partition function for this system may be written: 

[dU] e Ss 

[dU]Y,S A , e So (5) 

and lattice quantities may be evaluated as power series in 5 in the background S . This 
leads to a diagrammatic expansion related, but not identical to the conventional strong 



coupling j3 expansion [16| [[17]], the difference being that the strong coupling expectations 
are evaluated in a zero background. The actual diagrams used are also different, the first 
non-vanishing diagram in the strong coupling expansion for the mass gap being a closed 
cuboid of plaquettes, compared with the LDE, for which the first diagram is shown in 
Fig. 1(a). 

Calculation of the mass gap involves the evaluation of the connected correlation C(t) 
between two non-oriented plaquettes A and B with temporal separation ta in any spatial 
position. 

C(t) =< txU B (ta)tTU A (0) > c (6) 

The subscript C denotes the connected expectation or cumulant. The diagrammatic ex- 
pansion in powers of S has its first non- vanishing term at 0(5 l ). This is shown in Fig. 1(a), 
where a "ladder" of time-like plaquettes connects A and B. 

The next power in 5 adds one extra plaquette to Fig. 1(a) in all possible positions. Some 
examples of these are shown as Figs. 1(b)- (k). It should be noted that this calculation is 
carried out in the temporal gauge. This explains the absence of diagrams where a plaquette 
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is attached to the side of the ladder by a temporal link only. Such a link variable is set to 
unity, and therefore the extra plaquette is effectively disconnected. 

At this order, there is also a term proportional to S t So in the 5 expansion. This can 
be included as a partial derivative with respect to J of the 0(5 l ) diagram. 

Each of the diagrams shown has an associated multiplicity depending on its geometric 
properties. The basic diagram of Fig. 1(a) has a factor of 4Ri representing the fact that the 
ladder may be connected to any of the four sides of the lower plaquette A, which is taken 
as fixed, and that the upper plaquette B has R\ possible spatial orientations. Similarly 
the additional factor of 12R 2 in Fig. 1(b) arises from the R 2 possible spatial orientations 
of the extra plaquette, its possible attachment on any of the three sides of its neighbour, 
the fact that either of the two upper plaquettes could be B, and finally a factor of two to 
include the symmetrical configuration where the extra plaquette is attached to A instead. 
Note that Figs. 1(f), (h) and (k), which involve an additional plaquette in the body of the 
ladder, have a t-dependent multiplicity. 

Having enumerated the diagrams to the required order and calculated their associated 
multiplicities, their expectation values must be evaluated. The evaluation of simpler di- 



agrams consisting of up to four or five plaquettes by group integration |L3[ or character 



expansion |15[ has been discussed elsewhere. Another method is discussed in the appendix 
to the present paper. 

In the evaluation of the straight ladder diagram of Fig. 1(a) and subsequent modifica- 
tions thereof, an enormous simplification arises from the fact that the expectation value 
of a single link is a multiple of the identity (Eq. (Al)). This means that the contribution 
of Fig. 1(a) is just a product of factors representing the expectation values of the dou- 
bled links occurring on each rung. Remarkably this factorization extends to the connected 
expectation value, with the result that 

C 1(a) = (4Ri) 2V 2 6 V t+1 S t (7) 



Here each factor of V2 represents the expectation value of a single spacelike link of A and 
B, and we get a factor of V = dV2/dJ for each rung. The functions V n are defined as 
ratios of modified Bessel functions of argument 2 J: 

Two derived quantities which appear frequently in the contributions of higher order 
diagrams are V as defined above: 

V(J) = ^ (9) 

and 

Y(J) = AV A - 6K 3 + 2V 2 (10) 

The higher order diagrams consist of modifications to this basic diagram by inserting 
additional plaquettes at either end and/or in the middle. At order 5 t+1 only one of these 
alternatives is possible. The factorization property noted above extends to these higher 
order diagrams. That is, their connected expectation can be obtained from the basic 
building blocks shown in Fig. 2, with additional factors representing the bulk of the ladder. 

Diagrams which involve a modification at one end have a multiplicity which is inde- 
pendent of the total temporal separation t, whereas the three diagrams 1(f), 1(h) and l(k) 
which involve an addition to the middle of the ladder have t-dependent multiplicities. The 
latter essentially exponentiate in higher orders and so are the only ones which contribute 
to the mass gap when this is calculated from the ^-expansion of the ratio C(t + 1)/C(t) 
(See Eq. (|l6|)). The contributions of these diagrams are given below: 

Ci (f) = (12Rl(t - 1)) V t+2 V 2 8 5 t+1 

C m = {ARlit-l^V'V^Y -W 2 V)5 t+1 (11) 
C m = (2R x t) ^*-V 2 6 (Y 2 /3 - lQV 2 2 V 2 )5 t+1 
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The diagrams of order 5 t+2 are similarly built up by adding a further plaquette in all 
possible ways. Some examples are shown in Fig. 3. There are around 150 diagrams at 
this order, although again it is only those with t-dependent multiplicities which contribute 
to the ^-expansion of C{t + 1)/C{t). There are also additional terms in S l Sl and S t+1 S 
which arise from the expansion of the factor (S — So) r in Eq. ^. The most succinct way of 
including such contributions is to note that J always occurs in the combination J(l — 5) 
(See Eq. (||)). Thus the Vs occurring in the various expectation values are really functions 
of this argument, which needs to be Taylor expanded to the appropriate order. Altogether 
we may write 

C{t) = Y i M i D i (J{l-S)) (12) 

i 

where Mj denotes the multiplicity of the ith diagram and Di is its connected expectation 
value. 
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3. Extraction of the mass gap 

Having set up the diagrams necessary to calculate the correlation C(t), we now need to 
extract the mass gap m using the familiar result: 

C{t) ~ A e~ mat (13) 

as t — > oo, giving 

C(t) 

ma = lim In — (14) 

t-oo C(t + 1) V 7 

At first sight it might seem reasonable to calculate C(t) and C(t+1) separately, applying 
the PMS to each correlation, and then to extract ma from equation (0). However, this 
is not a fruitful procedure for two reasons. The first is that it is not in the spirit of the 
PMS, according to which it is the final quantity calculated which should be optimized 
with respect to J. More importantly, the convergence of the expansion performed in this 
manner is extremely slow. It is, after all, asking a great deal of a perturbation expansion, 
even when optimized, to give the correct t — > oo limit of C(t) with only a few terms of the 
expansion. 

The most important aspect of the problem is that some of the diagrams have multi- 
plicities which grow with t, reflecting the fact that additional plaquettes can be attached 
in a large number of positions to the body of the ladder. Thus the larger the value of t, 
the higher the order of the perturbation expansion required before the factorial denomi- 
nators in Eq. fl5|) eventually control the convergence. However, these diagrams essentially 
exponentiate. For example the series of "bracket" diagrams starting with Fig. 1(h) and 
continuing with Fig. 3(h) has the form of an exponential series for large t. Consequently, 
when the series for the ratio is taken the t-dependence cancels, as we show in more detail 
below. Thus by considering the Taylor expansion for the ratio, the limit t —>■ oo does not 
pose such a threat to the convergence of the series. A similar procedure was adopted by 
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Minister [I7| in the application of the strong coupling expansion to the calculation of the 
mass gap. 

We therefore apply the PMS to the Taylor expansion of the ratio C(t+ 1)/C(t), up to 
third order in 5. Writing the series for C(t) and C(t + 1) as 

C{t) = ^(6 + M + M 2 + ---) 
C(t+1) = 5 m (a + ai5 + a 2 5 2 + ...), (15) 

the ratio has the expansion 

C(t + 1) a 0jr , aib -a bi 2 6i(a 6i - aib ) - b (a b 2 - a 2 b ) 3 

^T = ^ + — + H 4+ - (16) 

This formulation leads to a naive large t limit for the mass gap. In going from temporal 
separation t to t + 1, we add an extra plaquette to the ladder part of each diagram, which 
gives an overall extra factor ~/3V(J) to the correlation. Thus one might expect the mass 
gap to be equal to — \n(^(3V{ J)). In fact, for the lowest-order contribution we have 

a = /oA m 

bo = /oA* (17) 

where A = \(3V , so that indeed a /bo = \(3V . In higher orders, however, t-dependent 
multiplicities give rise to corrections to the naive result. 
The next-order coefficients have the form: 

ai = (f 1 + (t + l)g 1 )X t+1 

bi = (fi + tg 1 )X t (18) 
Then to second order in 6, Eq. ( |T6| ) is: 

-^r =<5|A+<5 T 1 (19) 
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Again this is independent of t, and means that in this form of the expansion t does not need 
to be taken asymptotically large. It is sufficient to take it large enough for the diagrams 
to settle down to a generic form. 
At 0{5 t+2 ), writing 

a 2 = (f 2 + (t + l)g 2 + (t+l) 2 h 2 )\ t+1 

b 2 = (f 2 + tg 2 + t 2 h 2 )\ t , (20) 
the 0(5 3 ) term in Eq. (|16|) is: 



&i(a &i - ai&o) - b ( a o b 2 - a 2 b ) X 2 

O = 72 [fo{h 2 + 92) - fi9i + t(2h 2 f - gj] (21) 

°o Jo 

This apparently has a t-dependence, but in fact the coefficient h 2 is precisely \g\/fo 
because it arises from exponentiation of the t-dependent graphs at order 5 t+1 . As empha- 
sized by Minister |Tj], the summation over the spatial positions of the upper plaquette B, 



which also serves to project out zero spatial momentum in the correlator, is vital to this 
exponentiation. 

Altogether, then, we have the t-independent result for the ratio to order <5 3 : 

i— -| — - 

fo fo 



^- + ^ + ^ 2 (fo(h 2 + g 2 )-f 1 g 1 ) 



(22) 



C(t) 

This expression for C(t+1) /C(t) is still a function of J. According to the PMS criterion, 
we are looking for stationary points in J. Typical curves of the J dependence are shown in 
Figs. 4 and 5 for 0(S t+1 ) and 0(S t+2 ) respectively. At 0(S t+l ) there is a single maximum, 
and at 0(5 t+2 ) the value of C(t + l)/C{t) at the maximum is remarkably close to this, 
even though the position in J is quite different. 

At this stage our result is expressed in terms of the inverse of the lattice spacing a, 
which we need to take to zero in order to make contact with the continuum limit. The 
physical value of the glueball mass must in this limit become a fixed number times the 
SU(2) lattice scale Al- In the weak coupling limit, to two-loop level, this is given by 
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J ex P(-^rJ ( 23 ) 

We look for the constant C m such that 

ma = Crak^a , (24) 

which is the value for which the graph of Eq. ( |23"D against (3 is tangential to that of ma 
calculated in the LDE. The graphs are shown in Fig. 6 to 0(5 t+l ) and Fig. 7 to 0(8 t+2 ). 
These show good agreement between the orders, the tangents occurring at (3 — 2.62 and 
(3 = 2.64 at 0(5 t+1 ) and 0(S t+2 ) respectively. These results then give for the mass gap: 

m = 184A L {0(5 t+1 )) 

m = 197A L {0(5 t+2 )) (25) 



Compared to the strong-coupling expansion |17|, which gives m = 193A^ (f3 ~ 2.3) 
at order j3 6 and m = 127Al (j3 ~ 1.9) at order /3 s , our results show better consistency 
between consecutive orders; moreover, the /3-values where the tangents occur are further 
into the weak-coupling region. In a series expansion of this kind it is difficult to quote a 
precise error, but based on the difference between our two results at 0(5 t+1 ) and 0(S t+2 ) 
one would estimate the error as not more than ±13. 

Our results can be compared directly with those of Berg and Billoire ||18|| , who quote 
m = (190 ±10) Ax,. A comparison with the more recent work of Michael and Perantonis [IS 



on a 32 4 lattice is less straightforward, since they quote their results in lattice units and 
cast some doubt on the validity of asymptotic scaling. Nonetheless, converting m = 197A^ 
to lattice units at (3 — 2.5 gives ma = 0.70, in excellent agreement with their results. At 
(3 = 2.7 it gives ma = 0.42, which is slightly higher than their central value, but still within 
the error bars. 
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4. Conclusions 

In this paper we have demonstrated that the linear delta expansion with the principle of 
minimal sensitivity is a viable technique for the calculation of the mass gap for a lattice 
gauge theory. We have shown how the lattice diagrams appearing in this type of calculation 
can be easily evaluated by a process of building up chains of plaquettes from a simple 'root' 
diagram, and that connected expectations of these are as simple to deal with. The gauge 
fixing procedure adopted reduces the number of contributing diagrams, and makes them 
easier to evaluate. 

As always, the PMS is an integral part of the calculation. The potential ambiguity 
arising from the occurrence of multiple PMS points is not serious in this case. It is clear 
by comparison with the lower-order calculation that it is the broad maximum at 0(5 t+2 ) 
which is the appropriate one, and it is very encouraging that the resulting value of m is so 
stable in going from one order to the next. 

This calculation has shown the relationship between the LDE and the strong coupling 
expansion. The diagrammatic expansion used is similar, but the actual evaluation of the 
diagrams is different, requiring alternative techniques. 

It has proved sufficient to work with the correlators of simple plaquette operators rather 
than the more complicated "fuzzy" operators which have been found necessary in Monte- 
Carlo calculations. The fundamental reason for this is that we are effectively working 
on an infinite lattice, so that large separations are no problem, whereas in Monte-Carlo 
calculations it is necessary to enhance the signal at finite separations. 

The present calculation could be extended in various ways. In increasing order of 
difficulty these are: 

(i) to calculate higher mass glueball states. With a simple plaquette operator the J PC = 
2 ++ state occurring in the E + representation of the cubic group is accessible by weighting 
the different orientations of the upper plaquette. Other spin-parities would require larger 

13 



Wilson loops. 

(ii) to work with the gauge group SU(3) rather than SU(2). This would involve an extension 
of the techniques of the Appendix to SU(3). 

(hi) to go to next order in the 5 expansion. The difficulty here is the greatly increased 
number of diagrams which have to be taken into account and the consequent danger of 
missing an important contribution. 

Further possible extensions include calculations of the string tension and various quan- 
tities at finite temperature. Some work has already been done on these lines by Tan and 
Zheng p0 |, but using the free energy criterion mentioned above. It would be interesting 
to return to these problems using PMS optimization order by order in the quantity being 
calculated. 
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A Appendix 

In this appendix we wish to explain a method for calculating (connected) expectation 
values of a string of plaquette operators in the background of our trial action So (Eq. 
The method involves expressing the expectation values of products of single link operators 



in terms of tensor projection operators ||21|| , which can then be multiplied together within 
an algebraic manipulation package. Of those available, we found FORM the most suitable 
because it has an explicit summation convention and extensive substitution facilities. 

Let us start with the simplest example, < U >, where U is a single-link element of 
SU{2). It is clear that < U l a >oc (Y n ) l a ■= 5 l a , and by taking the trace we establish the 
coefficient as V 2 (see Eq. (||)): 



< lP a >= (V 2 Y n y a (Al) 

The notation has been chosen with a view to subsequent examples and in general refers to 
the Young tableau associated with a particular permutation symmetry, in this case trivial, 
of the upper indices. 

For the product of two C/'s belonging to the same link, there is similarly no difficulty 
(again by taking appropriate contractions) in showing that 

< UlUl >= (V 3 Y m + Y Q )Z (A2) 

where 

{Y^)l = \m + 5i5D 

(Y B y> b = \m - 5151) (A3) 

Here Ym and Yg are indeed projection operators which correspond to the two irreducible 
representations of the permutation group of the upper indices relative to the lower ones. 
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If we now go on to the product of three C/'s, there are three irreducible representations 
of the relevant permutation group S3. However, since we are constructing irreducible 
tensors of SU(2), the completely antisymmetric operator Yg effectively vanishes and only 
the completely symmetric and mixed symmetry operators Yj-m and Ym survive. 

The form of these tensors can be deduced from the group algebra of the conjugacy 
classes of S3, which comprise E, B := 3C 2 and C := 2C 3 . The non-trivial products in the 
algebra are B 2 = 3(E + C), BC = CB = 2B and C 2 = 2E + C. The vanishing of Yn 
for SU (2) corresponds to the equivalence relation B = E + C when acting on the identity 
element 5 l a 5 3 b 5 k . Thus in our search for projection operators we can limit ourselves to the 
sub-algebra of A 3 generated by the even permutations of E, C. It is then easy to construct 
the required projection operators as 

Y- = \{E + C) 
% = \{2E-C) (A4) 

acting on S l a 5l5^ by permutation of the upper indices. By taking traces we can establish 
the coefficients as V 4 , V 2 in the expansion 

< W a U{U k c >= (V A Y^ + V 2 Y W ) 5\5i5 k c (A5) 

For the product of four C/'s, we can anticipate that 

< UiuiU k c U l d >= (VsYcc- + VsY^n + Y m )8i8i8% (A6) 

The problem is to establish the specific form of the three projection operators. Again 
the procedure is to look at the group algebra of S4, which has five conjugacy classes. 
Because of the equivalence relations arising from the vanishing in SU{2) of completely 
antisymmetric combinations involving more than two indices, we can eliminate the classes 
of odd permutations and work with the group algebra of the alternating group A4. This 
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has three classes, E, B := 3C 2 , C := 8C 3 , with algebra B 2 = 3E + 2B, BC = CB = 3B 
and C 2 = 8E + 8-B + AC. From these it is possible to construct the projection operators 

V™ = -^(E + B + C) 

Y m = L(2(E + B)-C) (A7) 

These formulae are sufficient to evaluate all the diagrams we encountered up to 0(5 t+2 ). 
Diagrams involving W can be dealt with |yj by converting to U according to 

In particular, 

< Ufil j >= (\{V Z + l)y m + \{W, - 1)Y Q )% (A9) 

A given diagram will consist of a number of plaquettes with certain links in common. The 
procedure is then to write down the general expression for the corresponding amplitude, 
identify the shared links, apply the appropriate substitutions for < UU >, < UW >, 
< UUU > and < UUUU > and then sum over all repeated indices. 

In fact we actually want the connected expectation values (cumulants) < >c- These 
can be obtained by identifying the different terms in the expansion of Y[f =1 (d / dxi)h\Z 
with products of expectation values of the corresponding partitions of the iV plaquettes. 
Thus 

d d 1 d 2 Z 1 dZ 1 dZ 

corresponds to 

< UV > c =< UV > - < U >< V > (All) 

etc. 

These identifications can all be implemented within a short FORM program. 
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Figure 1: Non- vanishing diagrams at 0(5 l ) and 0(S t+l ) in the linear delta expansion for 
the 577(2) mass gap. The vertical direction represents time. The expression next to each 
diagram is its geometric multiplicity factor. The Ri are defined by Ri = 2d — 3 — i for d 
spacetime dimensions. 
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(b) 
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v 2 V 



V/(Y-8V 2 V) 



(c) 



V 2 (V 3 3 Y-8V 2 7 V) 
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V 2 2 (Y 2 /3-16V 2 'V) 



2 T t2 n 



Figure 2: Connected expectations of some of the basic diagrams used in calculating the 
correlation C(t). The symbols used are defined in the text. 
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Figure 3: Examples of non- vanishing diagrams at 0(5 t+2 ) in the delta expansion for the 
577(2) mass gap. The vertical direction represents time. Diagrams (a)-(c) have multiplicity 
independent of t, (d)-(f) have multiplicity oc t, and (g)-(i) have multiplicity oc t 2 . 



22 



PMS Graph at Order t+1 




Figure 4: 0(5 t+1 ) PMS graph {(3 = 2.62). 



PMS Graph at Order t+2 




Figure 5: 0(5 t+2 ) PMS graph {(3 = 2.64). 



Results to Order t+1 
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Figure 6: 0(5 t+1 ) results for ma (dotted line), compared with weak coupling renormaliza- 
tion group formula for A L a (solid line), with C m = 184. 



Results to Order t+2 
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Figure 7: 0(5 t+2 ) results for ma (dotted line), compared with weak coupling renormaliza- 
tion group formula for A^a (solid line), with C rn = 197. 
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